Particle beams guided by electromagnetic vortices: New solutions of the Lorentz, 
Schrodinger, Klein-Gordon, and Dirac equations 
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It is shown that electromagnetic vortices can act as beam guides for charged particles. The 
confinement in the transverse directions is due to the rotation of the electric and magnetic fields 
around the vortex line. Large class of exact solutions describing various types of relativistic beams 
formed by an electromagnetic wave with a simple vortex line is found both in the classical and in 
the quantum case. In the second case, the motion in the transverse direction is fully quantized. 
Particle trajectories trapped by a vortex are very similar to those in a helical undulator. 
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Electromagnetic waves with vortices have been ex- 
tensively studied both theoretically and experimentally. 
This field of research has became known as singular op- 
tics 0] . In this work I take these studies one step further 
and analyze the motion of charged particles in the vicin- 
ity of a vortex line. I shall consider the simplest possible 
solution of Maxwell equations with a straight vortex line 
and show that this configuration of the electromagnetic 
field acts as a perfect beam guide for charged particles. 
I study these nonspreading beams in the classical case, 
when the relativistic particle trajectory is determined by 
the Lorentz equations, and also in the quantum case, 
when the wave function describing the beam obeys the 
Schrodinger, Klein-Gordon, or the Dirac equation. In the 
classical and in the quantum case, I exhibit analytic so- 
lutions that enable one to fully understand the intricate 
dynamics of these beams. 

The electric and magnetic field vectors of my model 
Maxwell field are 

E(x, y, z, t) = Bquj (f(x, y, z, t),g(x, y, z, t), 0) , (la) 
B(x, y, z, t) = B k {-g[x, y, z, t), f(x, y, z, t), 0) , (lb) 

where B is the field amplitude measured in units of the 
magnetic field and 

f(x,y,z,t) = xcos(wt-) + ysm(ivt), (2a) 
g(x,y,z,t) — xsm(u)t-) — ycos(cL>i_), (2b) 

where t_ = t — z/c. This configuration of the field is the 
simplest example of the EM field with a vortex line (HQ. 
The solution of the Maxwell equations given by Eqs. (JTJ 
is not as artificial as it may look at a first glance. It 
is a fairly good approximation (near the z-axis and not 
far from the waist compared to the Raleigh range) to 
a realistic circularly polarized paraxial Laguerre-Gauss 
beam with n — and m = 1. 

The Lorentz equations of motion 



m("(r) = erK(r))6(r) l 



(3) 



for a particle moving in the field expressed in terms of 
the components (£, rj, £, 6) of the four- vector ^(r), have 



the form (for the sake of brevity, I shall occasionally drop 
the dependence on the proper time r) 



c6 



i = w c w/(f,»7,C,0)(0-C/c) , 
rj = ucug(t,ti,{,0)(0-t/c), 

c = —Um,v,c,o)+vg(t,v,c,o) 

c V 



(4a) 
(4b) 
(4c) 
(4d) 



where the dots denote derivatives with respect to r and 
lu c = eB /m is the cyclotron frequency. These equations 
are nonlinear but they can be explicitly solved owing to 
conservation laws. 

By subtracting Eq. (|4*c|l from Eq. (|4dl) . one obtains 
9 — C,/c = and this leads to the first conserved quantity 



6 - C/c = y l + (£ 2 + V 2 + C 2 )/c 2 -C/c = £ = constx. (5) 

Apart from the factor mc 2 , this constant is the light-front 
energy — the conjugate variable to t + = t + z/c 



l — v z /c y/ m 2 c 4 + p 2 c 2 — p z 



(6) 



Without any loss of generality one may assume that 
8(0) = = C(0) and then Eq. JSJ integrated with re- 
spect to t yields 



C/c = £t. 



(7) 



Thus, in this case, the proper time is proportional to the 
light-front variable. The second constant of motion is 
obtained by combining Eqs. H4al4cf> and it reads 

C-^(e 2 +?f) = ^-£)=const 2 . (8) 

Since the phase of the wave field changes in proper time 
with frequency io£ , I shall incorporate £ into tu and define 
the effective frequency £1 = £ui. 
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Owing to Eq. (J5J), the transverse motion separates from 
the longitudinal motion and the equations for £ and r\ 
may be solved first. This task is made easier by trans- 
forming Eqs. H4al4b|) to the frame rotating (in proper 
time) with the angular velocity f2/2 around the z-axis 
which amounts to replacing £ and r\ by the new variables 

a(r) = £(r)cos(ftr/2) +7?(T)sin(Qr/2) J (9a) 
P( T ) = -£(r) sin(ftT/2) +r?(r)cos(nr/2). (9b) 

The equations of motion for a and (3 read 

& = Q/3+ {Q 2 /4)a + uj c na, (10a) 
/? = -tta + (n 2 /A)(3- u c np. (10b) 

These equations result from the following Hamiltonian 

Despite the quadratic form of the Hamiltonian, it is 
not exactly a harmonic oscillator — the frequencies of 
the oscillations depend through £ on the initial condi- 
tions. Still, this Hamiltonian can be expressed in terms of 
the complex eigenmode amplitudes a± and a± (classical 
counterparts of the annihilation and creation operators), 

H = a*^a + — 17_ a*_a-, (12) 

where 57 ± = f2vl ± n/2 and k — Auj c /il is a dimension- 
less parameter that controls the particle behavior in the 
xy-plane. The amplitudes a± have the form 

a + = -./ — [pp-ja-i , (13a) 

2 V 7 V K + / 

a = l^(El±l^ +i(pa - J0) y (13 b) 

where 7 = \eBo\ and k± = vl±K. The minus sign in 
the diagonal form of the Hamiltonian l|12|) indicates that 
the beam dynamics in the transverse plane is governed 
by the same combination of the attractive/repulsive os- 
cillator forces and the Coriolis force as one encounters for 
a particle in the Paul trap |^|, an electron Trojan wave 
packet in an atom (cf. Eq. (7) of Ref. §), or Trojan 
asteroids in the Sun- Jupiter system ||. 

The general solution for £(t) and rj (r) is obtained by 
solving Eqs. I jlOjl in terms of eigenmodes and then undo- 
ing the rotation 10 ■ The final expression for the motion 
of particles in the plane perpendicular to the vortex line 
can be compactly written in the complex form 

£(r) + iri(r) = e inT ' 2 {{C + iDn+) sin(fi + r) 
-{Bk- - iA) sin(Q_r) - (D + iCk+) cos(0 + t) 

+ (Ak-- iB) cos(0_t)) , (14) 




FIG. 1: Two trajectories of electrons injected into the wave 
field Q with Bo = 10 -3 T and lo = 2tt ■ 10 9 s _1 . The initial 
longitudinal momentum p z of the electron is in both cases 
25keV/c but they have different transverse momenta. The 
narrow trajectory has p x = 5keV/c while the wide one (wig- 
gles on this trajectory are real) has p x — 50keV/c. The size 
of the box measured in wavelengths 2-kc/uj is 1^ x 1^ x 2|. 

where the constants A, B, C, D depend on the initial val- 
ues of the transverse positions and velocities 

A = (£0 + t?o/2cj c )/k_ , B = io/2w c , (15a) 
C = (770 + £o/2w c )/re+ , D = rj /2w c . (15b) 

For |«| < 1 one obtains bounded oscillations around the 
vortex line with four characteristic frequencies: f2 + ±f2/2 
and f2_ ± 51/2 and for |k| > 1 one has runaway solutions 
with exponential growth. The motion along the z-axis is 
obtained from Eq. JHJ by a straightforward integration. 
The resulting formula for £(r) has two parts — a part 
with oscillating terms and a linear part in r , 

C(r) = 

k 2 ui(D 2 -C 2 ) sin(2Q + r) k 2 lo(B 2 -A 2 ) sin(2fi_r) 

16ck+ 16ck_ 
_ k 2 lj CD(l-cos(2n+T)) k 2 uj AB(l-cos(2S7_r)) 
8ck + 8CK- 

+c i(\- £ + £j ^- {A2+B2+c2+D2) )- (16) 

Depending on the sign of the linear term, the guiding 
center of the beam may follow the electromagnetic wave 
or move in the opposite direction. By a special choice of 
initial conditions one may even get rid of the linear term 
altogether in which case the longitudinal motion will also 
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be bounded, but it requires fine tuning. This complex be- 
havior is a purely relativistic effect. In the nonrelativistic 
limit, the motion in the ^-direction is free, not affected by 
the wave at all, ((t) = v z t. In Fig.^ I show two trajecto- 
ries of electrons for different initial conditions. These tra- 
jectories are very similar to those in a helical undulator 
(an arrangement of permanent magnets used to produce 
circularly polarized radiation). In the present case, the 
role of permanent magnets is played by an electromag- 
netic wave with a vortex line and the beam confinement 
is due to a totally different (Trojan) mechanism. 

I shall start the analysis of the quantum-mechanical 
problem with the Klein-Gordon (KG) equation. The EM 
field may be derived from the vector potential 

A(x, y, z, t) = B (g(x, y, z, t), f(x, y, z, t), 0) . (17) 

As seen from the analysis of the classical solutions, it is 
preferable to use the coordinates x,y and t±. The KG 
equation in these coordinates reads 

where d± = d/dt±. Since the variable t + does not appear 
in this equation, one may seek its solutions in the form 

Mx,y,t.,t + ) =e-^ Mt ^ 2t -l M ^i>{x,y^). (19) 

An additional phase factor, dependent on i_, has been 
introduced to remove the mass term. The essential de- 
pendence on t— is still contained in the wave function ip. 
The function tp(x, y, i_) obeys the following equation 

+ ihfl c ((xd y + yd x ) cos(wi_)) ip 

- ihil c ((xd x - yd y )sm(ujt -)))$, (20) 

where Aj_ is the transverse part of the Laplacian and 
Cl c = eB /M. This equation is exactly the same as a 
nonrelativistic Schrodinger equation except that the role 
of the mass m is played by the separation constant M and 
the time parameter is replaced by the light-front variable 
t_ . Therefore, everything that one can say about the 
solutions of the Eq. 1)20(1 applies to the solutions of the 
Schrodinger equation. Upon transforming the equation 
(I20|l to a comoving frame by the substitution 

ip = exp(-^-(xd y - yd x ))<p, (21) 
one finally obtains 

+ih((n c + u/2)xd y +(n c -u/2)yd x )(f>. (22) 



By rearranging the terms, one may establish that the 
particle in this frame moves effectively under the influ- 
ence of the constant magnetic field B = (0,0,mw/e) 
and an additional repulsive quadratic potential V = 
-(Muj 2 /8)(k 2 + x 2 + nly 2 ). 

All stationary solutions of Eq. (1221) are most easily clas- 
sified with the use of the creation and annihilation oper- 
ators. These operators diagonalize the Hamiltonian 

H = {pi +%)/2M + Mn 2 (x 2 + y 2 )/2 

-(O c +uj/2)xp y - (O c - uj/2)yp x (23) 

and are obtained from the classical amplitudes ((13(1 by 
the replacements 

a± — > Vha± , Oj. — > y/ha/j- , (24a) 
{at,P,p a ,Pp)-> (x,y,p x ,py). (24b) 

This leads to the following form of the Hamiltonian 

H = huj (n + (a* + a+ + 1/2) - K + (a*_a- + 1/2)) . (25) 

Thus, in contrast to the Volkov solution in the plane 
wave EM field 0, the motion in the transverse direc- 
tion is fully quantized. In contrast to the motion in a 
constant magnetic field, the particle is localized near the 
z-axis. Different normalization of the classical and quan- 
tum Hamiltonian is due to the fact that the first one gen- 
erates the evolution in proper time, while the second one 
generates the evolution in the i_ variable. These two pa- 
rameters differ by the scaling factor £ . The quantum the- 
ory becomes consistent with the classical one when M/m 
is identified with £. It means that M/m = Q/uj and, as 
a result, the value of k encountered in quantum theory 
becomes equal to the classical one (Aui c /Q, = 4J7 c /oj), as 
it should be. 

Having diagonalized the Hamiltonian, one may gener- 
ate the whole Fock space of stationary solutions. They 
are obtained by acting on the fundamental state <po with 
the creation operators. The fundamental state is the one 
annihilated by both operators a±. Solving two simple 
differential equations a±cf>o — 0, one obtains 

M*,V) =NeM~x 2 /d 2 + ~ y 2 /dl - ixy/d 2 ), (26) 
where the parameters d± and d are given by 
= hl + K+ K- d2 = fi ! + «+«- 

7 K± 7 1 — 

The wave functions of the Fock states are polynomials in 
x and y multiplied by the Gaussian 1(261) . In the labo- 
ratory frame these solutions are not stationary since the 
beams do not exhibit rotational symmetry around the z- 
axis (d+ =/= d-). In particular, the fundamental solution 
takes on a form of a rotating helix. 

There is also a plethora of nonstationary solutions of 
Eq. ((22(1 . First, there are those that correspond directly 
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to classical trajectories — the analogs of coherent states. 
The fundamental solution 1|26|) corresponds to a trajec- 
tory which just sits on the vortex line, but one may easily 
obtain solutions of the KG equation representing all other 
classical trajectories. According to a general scheme Q 
valid for all quadratic Hamiltonians, displacing any solu- 
tion of the KG equation by the solutions of the classical 
equations of motion leads to new solutions. Applying 
such displacements to the solution one obtains 

0(x,V,t-) = N{T)e l{xPa{T)+vp ^ T))/h 

x <j> (x-a(T),y-P(r)), (28) 

where r = t_/£ and the center-of-mass trajectories are 
obtained by solving the Hamilton's equations of motion 
that follow from (fill) . The time-dependent phase of the 
normalization constant is equal to the classical action ■ 
To obtain the solution of the original equation, one must 
transform the wave function from the comoving frame 
back to the laboratory frame applying the inverse trans- 
formation to i|21|) . Only then one obtains the quantum- 
mechanical counterparts of the classical trajectories. 

Solutions based on a rigid Gaussian — the analogs 
of coherent states — do not exhaust all possibilities. 
Since the center of mass motion of the Gaussian wave- 
packet decouples from its internal motion, one may easily 
generate solutions based on oscillating Gaussians — the 
analogs of squeezed states. The Gaussian parameters d± 
and d for such states are functions of t^. These states do 
not have direct classical counterparts and their complete 
analysis will be given elsewhere. 

The solution of the Dirac equation proceeds along sim- 
ilar lines. I begin with rewriting the Dirac equation in 
the electromagnetic field (1) 

ihOt^ = (ca-(-ihV - eA) + (3mc 2 ) (29) 

as a set of two coupled equations for the two-component 
wave functions 

2ihd + ^ + = c(mco- z - <r^-(ihV + eA)) (30a) 
2ihd-^- = c (mca z - er ± • (ihV + eA)) * + , (30b) 

obtained with the use of the projections P± — (1 ± a z )/2 

* = (P+ +P- )* = *+ + (31) 

The dependence on the variable t + can again be sepa- 
rated by the same substitution lj2*T)l . leading to 

Mc& + =(mca z -<Tx-(iKV + eA))V-, (32) 

(Til/ 'Cp 1 \ ~ 
2ihd- +— — J = c(mc(T z - crj_-(ihV + eA))W+, 

(33) 



The first equation enables one to express \E'+ in terms of 
ty- and leads to Eq. (|20|l for Again, as in the case 
of the KG equation, the dependence of the potential on 
t- may be eliminated by the substitution l|19f) and the 
equation for can be reduced to the same Eq. (|2*2l as 
for a spinless particle. Still, the spin does play a role in 
the Dirac particle dynamics. Since the transformation to 
the comoving frame should also involve the spin part, the 
proper transformation rule, instead of 1)21(1. is now 

*± = cxp{-^-(xd y - yd x + ia z /2))$±. (34) 
Finally, the wave equation for in the comoving frame 



ihd-<f>- = 




+ih ((n c + ^)xd y + (n c - f (35) 

differs from Eq. 122|) only by a simple spin term. Every- 
thing that has been said before about stationary solutions 
of the KG equation applies with almost trivial changes 
to the Dirac equation. 

There are two properties of the solutions of the wave 
equations described here that might lead to new effects: 
the quantization of the transverse motion and the break- 
ing of the rotational symmetry. This may help to observe 
the rotational frequency shift predicted some time ago || 
that depends crucially on these features. 
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